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Abstract 

Using the coupled cluster expansion with the random phase approximation, we calculate the 
long wavelength vacuum wave function and the vacuum energy of 2+1 dimensional Hamiltonian 
SU(2) lattice gauge theory (LGT) up to the seventh order. The coefficients fiQ , /i2 of the vacuum 
wave function show good scaling behavior and convergence in high order calculations. 
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1 Introduction 



Lattice gauge theory was developed into a promising first principles approach to the nonpertur- 
bative aspects of gauge field systems. Most of our present knowledge about LGT has obtained 
from numerical simulations. However, in order to gain more physical insight into the theory, it is 
desirable to develop more analytical methods. Based on Greensitc's proposal [1], we developed the 
coupled cluster expansions for the Hamiltonian LGT [2, 3, 4]. Although the preliminary research 
shows that this scheme is reasonable to study the LGT, it suffers from the rapid proliferation and 
non-independent of clusters produced in the high order expansive calculations. 

In Ref.[2], we used the Cayley-Hamilton relations 



to eliminate redundancies and find the independent cluster bases, where U is any group element. 
But those relations are too complicated to be used in the high order expansions. We have tried to 
use the improved Hamiltonian [5] with the tadpole improvement [6], so that the physical results 
could be obtained in a relative low order expansion. The research results for 2+1 dimensional U(l) 
LGT shows that the improvement of convergence is immaterial comparing with the complex of 
calculation brought about by the improved Hamiltonian [7] . Recently, we introduced the random 
phase approximation (RPA) into the coupled cluster expansions to circumvent the above problems, 
the preliminary investigation results were encouraging [8]. In this paper, we use this new method 
to calculate the vacuum wave function and vacuum energy of 2+1 dimensional SU(2) LGT. 

The paper is organized as follows. In Sec. II, we briefly review the scheme of the coupled 
cluster expansion with RPA. Section III is devoted to the calculation of the vacuum wave function. 
In Sec. IV, the conclusions and discussions are presented. 

2 Formulation and Approximation 

The Kogut-Susskind Hamiltonian [9] is 



TrU = TrU+, = 



UTrU-1 for SU(2) 



(1) 



2TrU+ = {TrU f - TrU^, = U^TrU - UTrU+ + 1 for SU(3) 



(2) 




(3) 
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where the index I denotes the hnks between sites, a is the lattice spacing, and Up the plaquette. 
The vacuum wave function can be written as [1] 

V'o(C/) = e^(^). (4) 

where R{U) consists of various Wilson loops or linked clusters with appropriate symmetries for 
the state. In the continuum limit(a 0, or equivalently P = —>■ oo), the long wavelength 
behavior of the vacuum state can be approximated by [2] 

ip{A) = exp (^-fio J (fxtvF'^{x) - 1^2 J (fxti[DiF'^{x)] + • • . (5) 

with F being the field strength tensor and Di the covariant derivative. The superrenormalizability 

of the theory in 2+1 dimensions implies that fiQ — > const. /e^ and ^2 const. /e^, where e is the 
invariant charge which is related to the dimensionless coupling constant g hy = e^a. 

The eigenvalue equation for the vacuum state of H is 

[E^, RiU)]] + [E^, R{U)m, R{U)]} - ^Y.^iUp = ^eo, (6) 

I y p y 

where eo is the vacuum energy. In the coupled cluster expansion, R is expanded in terms of a set 
of linearly independent cluster Gn,i with suitable symmetries 

M 

R{U) = Y,Rn = Y.Cn,iGn,i, (7) 
n=l n,i 

with Gn,i denoting the i-th cluster of order n, c„ j being a coefficient to be determined by Eq.(6), 
and M the highest order number in the coupled cluster expansion. Substituting it into Eq.(6) and 
adopting the truncation scheme of Ref. [2, 3, 4], we obtain the truncated eigenvalue equation 

M 



r 14 2 

Y,\[Eu[EuY,Rnm]+ E [EuRn{U)][Ei,Rr,\U)\\--Y^iiUp = —eo. 

I ^ n=l n+n'<M ^ 9 p 9 



(8) 



According to the symmetry of the vacuum state, the lowest order term Ri is chosen to be 
composed of just one cluster, which is an elementary plaquette 

i?i = ci,iGi,i = ci,iEtrC/p. (9) 
p 

The term [Ef,Rn][Ef,Rn'] in Eq.(6) will produce the (n + n')-th order clusters. For example, 
[E'l , Ri][Ef , Ri] generates the second order clusters involving two plaquettes. Obviously, clusters 
with order (n + n') may contain at most {n + n') Wilson loops. But not all clusters produced in 



such a way are independent. They may be related to each other by relation (1). We have to use 
this relation to identify and get the independent clusters. When clusters have many Wilson loops, 
relation (1) becomes so complicated that it is impossible to find a set of independent clusters, as 
pointed in Sec. I. To avoid this difficult, the RPA is applied in the expansions [8]. In RPA, one 
set of operator in a product of two sets of operators is replaced by its average value [10]. Here, we 
replace one Wilson loop by its vacuum average value when a cluster produced by [Ef, Rn][Ef, 
consists of two Wilson loops. It is easy to prove that any cluster produced in the above procedure 
contains only one Wilson loop. Suppose Rn and are linear combinations of clusters which 
contain only one Wilson loop, then, the new clusters produced by [Ei, Rn][Ei, Rn'] will contain 
at most two loops, and one loop will be replaced by its vacuum average after applying the RPA. 
On the other hand, according to Eq. (9) the first order cluster Gi^i has only one Wilson loop. 
Therefore, all new clusters produced in the calculation consist of only one Wilson loop after using 
the RPA. 

Because the clusters consist of only one Wilson loop, the independent cluster bases for the 
expansions can be obtained directly. In addition, the number of independent bases of high order 

expansions is much smaller than that without using the random phase approximation, for example, 
the number of independent third order clusters is nine in the coupled cluster expansion [2], while 
it is two after using the RPA ( see Sec. Ill ) . Therefore, the calculation is simplified very much. 

The vacuum average value of a Wilson loop can be determined by the Feynman-Hellman the- 
orem. Let G be some Wilson loop and < G > be its vacuum average. Defining W = H2a/g^, we 
make the following change [11]: 



W 



(10) 



where is a variable and will take zero at last. Prom W'^\ipo >= Wq\iI>o >, we get 



<G>= 



(11) 



3 Calculation of approximation 



We now present the calculation of the expansion. Because 



[Ef, Gi,,][Ef, Gi,i] = -4 - 2G2,i + G[ + G'^ 



2) 



(12) 
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there are three new clusters produced in the second order calculation, where the three clusters 
denoted by ^2,1, G[, and G'2 are given in Fig. 1. Two one which consist of two Wilson loops turn 
to cluster Gi^i times < Gi^i > by RPA. Therefore, there is only one cluster at second order after 
applying the PRA, that is ^2,1, and 

R2 = C2,lG2,l. (13) 

Substituting Ri and R2 into Eq. (8), we obtain a set of equations about ci,i, C2,i and wq with 
a parameter < Gi,i >. The parameter < Gi^i > can be determined by Eq. (11). Solving those 
equations, we get the second order approximation of vacuum wave function ^^^{U) ss qRi{U)+Ri{U) _ 
The long wavelength coefficients up to the second order are 

M = [^ + 2c2,i]/, (14) 
= -'^9'. (15) 

The third order clusters are produced by term [Ef , Ri][Ef , R2], i.e. ci,iC2,i[£^f, G2,i]. 
Since 

1 3 

[Ef, Gi^iWEf, G2,i\ = -SGsj - G3,2 + ^3 + -G4 + -G5 - 3Gi,i 

~ -2G3,i - G3,2 + 3 < Gi,i > G2,i - 3Gi,i, (16) 

we get 

-R3 = C3,lG3,i + C3,2G3,2. (17) 

There are only two independent clusters at order 3, while it is nine in the expansion without RPA 
[2], so the calculation is much simpler. In Eq. (16), when applying RPA to cluster G3, G4 and G5, 
we replace the small Wilson loop with its vacuum average and let the large one remain unchanged 
as in Ref. [8]. Thus, only one vacuum average of cluster emerges in the third order calculation. 
This makes the calculation simple further. On the other hand, the vacuum or exciting states 
possess definite correlation lengths. Only when the space occupied by glueball is covered with the 
Wilson loops, the calculation is efficient, so, we replace the small loop with its vacuum average 
and preserve the large loop. Submitting R2 and i?3 into Eq. 8, and solving the equation, we 
obtain the third order approximation of the vacuum wave function. 

Higher order calculation can be carried out similarly. We have done the calculation up to the 
seventh order, the number of independent clusters are 1, 1, 2, 6, 14, 44, and 109 at order 1, 2, 3, 
4, 5, 6, and 7 respectively. The results for ^0 and 1x2 from the fourth order to the seventh order 
are presented in Fig. 2. In Fig. 3, the vacuum energy wq against are plotted. 



4 Results and discussions 



Prom Fig.2, we see that the curves of /xq (or /X2) show good scahng behavior and convergent trend 
in weak coupHng region /3=5.5 to 9.5. The sixth and seventh order values of //q are coincident in 

the scaUng region, which shows that the approach is effective and rapidly convergent. From the 
seventh order results, we obtain 

fio = 2.3, (18) 
H2 = -0.3. (19) 

A Monte Carlo measurement was gave by Arisue [12] 

Ho = (0.91 ±0.02), 
//2 = -(0.19 ±0.05). (20) 

The value of /xq in Eq. (18) is larger than Arisue's value. It is perhaps that there is some systematic 
error introduced by the random phase approximation in our procedure. But it is also probable 

that the Monte carlo results of Eq. (20) was not enough accurate. Because the values in Eq. (20) 
can be reproduced in the low order ( the third order ) expansion without RPA, while in higher 
order ( the fourth order ) calculation the values became large [2]. 

In Fig. 3, the curve of the sixth and seventh order results are almost coincident, which proves 
the expansion is able to converge rapidly again. We also give the third order result of the vacuum 
energy calculated without RPA for comparing. The third order values of the vacuum energy with 
RPA are lower than that without RPA. Such a case was also true in the case of SU(3) LGT [8]. 
We think it means that the vacuum state obtained by using RPA is better than that without using 
RPA, and it may explain the rapid convergence of the approach. 
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Figure Captions 

Fig.l The linked clusters used in the expansions up to order 3. 

Fig. 2 /iQ and ^2 as a function oi (3 = A/ . The four curves represents the results from the fourth 
order to the seventh order expansion with RPA respectively. 

Fig.3 The vacuum energy vs (3. We also give the third order vacuum energy calculated without 
using RPA. 
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Fig. 1 The linked clusters used in the expansions up to order 3. 
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